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Intensity interferometry of thermal photons, having transverse momenta kx ~ 0.1 - 2.0 GeV, 
produced in relativistic collision of heavy nuclei is studied. It is seen to provide an accurate infor- 
mation about the temporal and spatial structure of the interacting system. The source dimensions 
and their kx dependence revealed by the photon interferometry, display a richness not seen in pion 
interferometry. We attribute this to difference in the source functions, the fact that photons come 
out from every stage of the collision and from every point in the system, and the fact that the rate of 
production of photons is different for the quark-gluon plasma, which dominates the early hot stage, 
and the hadronic matter which populates the last phase of the collision dynamics. The usefulness 
of this procedure is demonstrated by an application to collision of lead nuclei at the CERN SPS. 
I/-} [ Prediction for the transverse momentum dependence of the sizes for SPS, RHIC, and LHC energies 

. are given. 
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I. INTRODUCTION 



The search for quark hadron phase transition and the investigation of the properties of quark gluon plasma stand 
among the most challenging as well as rewarding pursuits of high energy nuclear physics today. The observations 
of jet-quenching P, @ and the elliptic flow 0-13 ■ and the success of the partonic recombination as a model for 
hadronization in recent experiments at the Relativistic Heavy Ion Collider at Brookhaven presage the shift of the 
focus to even more interesting questions about how the plasma is formed and how does it evolve. 

It is expected that the quantum statistical interference between identical particles emitted from the relativistic heavy 
■^j- ion collisions will provide valuable inputs for these investigations. The pion intensity interferometry experiments have 
however thrown up some startling questions, which need to be answered before we can address these issues. The 
hydrodynamics models which have traditionally provided a quantitative description to data on particle spectra and 
elliptic flow, fail to explain the near identity of the so-called outward and side-ward radii seen in these experiments 
and invariably lead to outward radii which are up to a factor of 2 larger than the side- ward radii [|| . 

An investigation of the structure and dimensions of the source in terms of photon interferometry should provide 
a more direct and valuable insight into these questions. The advantages of using photons for such studies is well 
known 0; they interact only weakly with the system after their production and are free from such distorting effects 
as the re-scattering || and Coulomb interactions, which have dogged these investigations using identical hadrons. 
Photons are also emitted at every stage of the collision dynamics and from every point in the system- and not only 
from the freeze-out surface. We note though that some studies have postulated that pions may be emitted not only 
from the surface as normally assumed; they also escape the system continuously from all points in the system with 
some escape probability. A large probability for this will however be hard to justify, as the pions do interact strongly 
with the hadroinc medium. Still this possibility along with the modifications introduced due to final state interactions 
for hadrons, necessitates a great deal of further study before the radii obtained from the intensity interferometry of 
hadrons are used to get the spatial and temporal extent of the source in these collisions. 

Theoretical investigations of the nature of the photon correlation function for relativistic heavy collisions, where a 
quark-gluon plasma may be formed have been carried out by several authors [ToL ITU . An experimental measurement 
is a far more difficult proposition though, due to the huge back-ground of decay photons and a considerably smaller 
production of direct photons. Thus, so far only one experiment at cyclotron energies |l2j | ( where direct photons mostly 
originate from bremsstrahlung of protons against neutrons) , and a measurement by the WA98 collaboration [p| (where 
photons originate from hadronic reactions and possibly quark matter) for the central collision of lead nuclei at CERN 
SPS have been reported. 

It is expected that such experiments will have a higher possibility of success at RHIC and LHC energies. Firstly, 
the larger initial temperatures expected there will lead to a larger production of direct photons. Secondly, a large 
suppression of high momentum pions due to the onset of jet-quenching will lead to a reduction in the back-ground from 
the decay photons. Furthermore, a large production of photons having high kx is expected from the pre-equilibrium 
stage of the collision |Tj, [l5| , which can be reliabl y cr afted using the parton cascade model [16| and it has been shown 
to have a very distinct distribution of the source \u\ . 

The early theoretical investigations 0, |n| of photon interferometry were also aimed at isolating photons coming 
from the quark matter from those coming from the hadronic matter. These studies concentrated on photons having 
high transverse momenta (kr 3> 1 GeV), which are more likely to be emitted from the early hot and dense phase of 
the plasma. 
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In the present work we focus our attention on the photons having low and intermediate for ~ 100 MeV - 2 GeV. 
The photons at the lower end of the fcp range considered would have a dominant contribution from the late stages of 
the hadronic phase. This should reveal a source which is strongly affected by the radial flow. The higher end of the 
kx should have a leading contribution from the early hot and dense stage of the plasma, when the flow is still small 
and the temperature is large. 

We also incorporate several improvements over the early exploratory works [lfj . Firstly, we use the complete leading 
order results for the production of photons from the quark matter 18], which has important contributions from the 
bremsstrahlung, and annihilation of off-shell quarks [l9( ■ Secondly, for hadronic reactions we use the state of the art 
results from Turbide et al. |2(J] along with inclusion of hadronic form-factors at the vortexes, incorporating strange 
mesons. The dominant hadronic bremsstrahlung process (itit — * itttj) is included |2ll] for low fop photons and double 
counting via tttt — ► pry and p — > ttttj avoided by limiting the contributions of the latter two processes to E^ > 500 
MeV. This should be adequate till a more detailed and complete calculation is available. Finally, we use a much richer 
equation of state for the hadronic matter, with the inclusion of all the particles in the particle data book, having M < 
2.5 GeV. Results of our model calculations are compared with the recent data obtained by the WA98 experiment and 
predictions given for RHIC and LHC energies. 

We shall see that the competition between the cooling due to expansion and the blue-shift of the photon spectra 
due to transverse expansion give rise to a unique fcp dependence of the outward and side-ward correlation radii, which 
evolves rapidly as the initial temperature of the system changes, as we go from SPS to RHIC to LHC. 

We first discuss the tools necessary for the analysis of the intensity interferometry and briefly discuss the initial 
state and the rate of production of photons. Next we check whether the description used by us provides a reasonable 
description of the single photon spectra measured at . The results for some typical momenta of the photons and 
transverse momentum dependence of correlation radii are discussed next. Finally we give our conclusions. 



II. FORMULATION 

A. The correlation function 

The spin averaged intensity correlation between two photons with momenta ki and k 2 , emitted from a completely 
chaotic source, is given by: 

1 |J d 4 zS(x,K)e"-«| 2 
C(q, K ) - 1 + 2 J d 4 x S(x, kx) / d*x S(x, k 2 ) (1) 

where S(x, k) is the space-time emission function, and 

q = k!-k 2 , K = (ki+k 2 )/2. (2) 

In a hydrodynamics calculation, the space-time emission function S is replaced by the rate of production of photon 
EdN/d ,xd k, in the hadronic or the quark matter, as discussed earlier. 

The results for the correlation function C(q, K) are best discussed in terms of the so-called outward, side- ward, 
and longitudinal momentum differences. Thus for the four-momentum k^ of the zth photon, we have, 

fof = (hr cosh y t ,ki) (3) 

with 

ki = {k iT cos'ipi, k iT sin^i, for sinh yi) (4) 

where kx is the transverse momentum, y is the rapidity, and ip is the azimuthal angle. Defining the difference and 
the average of the transverse momenta, 

qT = kiT - k 2 T ,K T = (kix + k 2T )/2 , (5) 
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The corresponding radii are obtained by approximating: 



(6) 



(7) 



(8) 



C(?out, tfside, 9long) - 1 + jT CXp [- (^ ut -Rout + ^sido-^sidc + ftongAong) / 2 ] 



(9) 



Note that for this choice of parametrization for the correlation function J2| = 1/ < q\ >, where i = out, side and 
long, and the average is performed over the distribution (C — 1). 

A one dimensional analysis of the correlation function C is some-times performed in terms of the invariant momen- 
tum difference 



1 



C(JqL) = 1 + 7y cxp 



(10) 



where 



= \j2k 1T k 2 T [cosh(yi - j/ 2 ) - cos(V>i - ip 2 )]. 



(11) 



However the significance of the corresponding radius does not have a clear meaning |23j, and it is also not useful 
for comparing results obtained using different particles. The variable q 2 + q 2 = qf nv + 2q 2 is also occasionally used for 
a one dimensional analysis p4j|. and one writes: 



C^^ + q 2 ) = 1 + \ CX P [-(q 2 + q 2 )R 2 /2] 



(12) 



lout + <Zsidc + C J\m 



Note also that q 2 

Here, it is useful to realize that in actual studies of photon interferometry the correlation function C would be 
parametrized as: 



C(<2out, gsidc, <7long) — 1 + - X A X exp [- (<7 2 ut i? 2 ut + g 2 ido^ 2 idc + ftong-Rfong) / 2 ] 



(13) 



where A is a measure of single (S) vs. decay photons (D), under the assumption that the source is completely chaotic. 
Obviously the later do not contribute to the correlation (the life-time of it is ~ 10~ 16 seconds) and thus: 



C(<jout = 0, <7 S idc = 0, <7io„ g = 0) = 1 + -A 



so that 



A 



S 2 



(S + D) 



(14) 



(15) 



and therefore using the measured values of the total photon yield (S + D) at the momentum K, we can get the yield 
of single photons from the intercept of C on the y-axis (Peressounko 1 1 1| ^ - This can provide a check on the results on 
single photon obtained by, say, a subtraction of photons from decay of pions and etas. 
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FIG. 1: The rate of emission of single photons from hadronic reactions considered in Ref. [20l| at T — 180 MeV, with the 
inclusion (solid curve) of hadronic vortex form-factors. 
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FIG. 2: (Colour On-line) The production of thermal photons from the quark matter plus hadronic matter, with and without 
the inclusion of the hadronic form-factors in central collision of lead nuclei at CERN energies. 

It is important to realize that this procedure is not directly applicable to the case where only a one-dimensional 
analysis in terms of q- lnv (Ealll)[l is performed. The photon-pairs having transverse momenta in some bin around K^, 
rapidity difference within Ay, and azimuthal angle difference within Aip will lead to a range of outward, side-ward, 
and longitudinal momentum differences which will decide the actual correlation function C (Eq[!5J) . Even a given value 
of gi n v will admit a range of these momentum-differences, and thus the value of one-dimensional C will be obtained 
by taking an average over the corresponding results. This will lead to an effective A which will not reduce to unity, 
even if there were no decay-photons [23. This "effective" A, as well as the i?i nv will depend upon Kt, AKt, Ay, Aip, 
and the correlation radii, i? ut, ^sidc, and i?i on g- Thus the direct photon yield can be measured only by assigning 
a value to these radii and becomes model dependent. Of-course, i?i nv can be assigned, though its meaning remains 
unclear and it remains model dependent. 

A simple instance should exemplify this issue. Consider a case where 2/1 = 2/2 = and ip± = ip2 so that g s ide = 
gi ong =0 (see Eq. H7l8|l ). In this case q out = (/cit — k 2 r), but q- lm = 0. Thus the "true" correlation function (Eq^J) 
will be decided by the value of the outward momentum difference and the outward correlation radius, while naively 
one would expect the one dimensional correlation function (Ea 110(1 . to reduce to 1.5. Thus even for a case when there 
are no decay photons, we would need a correction factor to account for the range of values admitted by the "true" 
correlation function. 

The procedure may however still work for the one-dimensional analysis in terms of the variable \/(/q + q 2 (Ea ll2fl . 
as this vanishes only when all the three momentum-differences, <7out, 9sidc and q\ on g also vanish simultaneously. The 
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FIG. 3: (Colour On-line) The production of single photons in collision of lead nuclei at CERN . 

meaning of the corresponding radius and its relation to the outward, side-ward and longitudinal correlation radii, 
remains unclear though. 

While we are discussing the merits of the one dimensional analysis in terms of the invariant momentum differ- 
ence, gw, we recall that for typical cases, the correlation function C differs substantially from unity for momentum 
differences < 0.2 GeV. For studies using single photons this limit has to be less than m^o, to stay clear of the 7r° 
peak. If kiT are large, choosing Ay and Aip as small, we can get very small values of the side-ward, outward, and 
longitudinal momentum differences. However as gw oc kx, small values for gj nv can be obtained only if kiT are quite 
small. For larger value of kiT, the necessary Ay and Aip bins would be too small to admit meaningful statistics. 
Thus, it is no wonder that both the photon interferometry experiments reported in the literature so far [3 13, use 
the one-dimensional analysis in terms of gi nv and utilize photons having very low /ct- We hope that a much larger 
statistics expected at RHIC and LHC will help us get over this problem by providing results for the full 3-dimensional 
correlation function. 



B. Single photons 

We very briefly recall the treatment for the space-time evolution of the system and the production mechanism of 
the photons in these calculations. 

For the first set of calculations, we consider central collision of lead nuclei, corresponding to the conditions realized 
at the CERN SPS. We further assume that a thermally and chemically equilibrated quark-gluon plasma is produced 
in such collisions at the initial time r , and use the assumption of is-entropic expansion to estimate the initial 
temperature, To. Thus we have 

where At is the transverse area of the system, dN/dy is the particle rapidity density, and a = 42.257r 2 /90 for a 
plasma of massless quarks (u, d, and s) and gluons. The number of flavours for this purpose is taken as w 2.5 to 
account for the mass of the strange quarks. We further assume a rapid thermalization of the plasma, constrained by 
the uncertainty relation, r = 1/3T . The initial energy density is taken as proportional to the "wounded - nucleon" 
distribution, appropriate for SPS energies. 

The quark-hadron phase transition is assumed to take place at 180 MeV, and the freeze-out at 100 MeV. The 
relevant hydrodynamic equations are solved under the assumption of boost-invariant longitudinal and azimuthally 
symmetric transverse expansion using the procedure discussed earlier [2f| and integration perf ormed over the history 
of evolution. This procedure is known to give an accurate description of hadronic spectra |27| . 

We have already indicated the sources for evaluation of rates of photon production from quark and hadronic matter. 

As a first step we plot the rate of production of single photons due to all the reactions among mesons considered 
in Ref.|2£j. (see Fig. ^| at T = 180 MeV. We see that the inclusion of the form factor considerably reduces the 
rate for production of photons having higher energies, leaving the results for the lowest energy photons unaltered. 
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We also give the results for the rates from quark matter at the same temperature. It should be remembered that 
photon intensity interfcrometry should be less sensitive to the specific details of the rates of production from quark 
and hadronic matter, but should depend on their relative strengths. 

In order to uniquely establish the importance of the mandatory requirement of including the form-factors discussed 
in Ref. |2GJ, we plot (see Fig- the sum of the production of thermal photons from the quark matter and the hadronic 
matter in a collision of lead nuclei at the CERN SPS, with the initial conditions discussed earlier. The data obtained 
by the WA98 experiments 0, and the production from the quark-matter are also shown for a comparison. We 
immediately note that the exclusion of the form- factors will considerably overestimate the production of single photons 
at higher kx, as it is well established from several studies that up to half of the single photons measured by the WA98 
experiment have their origin in the prompt-QCD process, which should be added to the thermal photons to get a 
quantitative description of the production of the single photons. We further note that the calculations underestimate 
the data by a factor of 7 - 9 at the lowest kx reported recently (see discussion later). 

Finally, in Fig. [3] we show our results for the production of single photons, by summing the thermal contributions 
(with the inclusion of form-factors while estimatin g th e hadronic rates) and the prompt-QCD contribution. For the 
later we use a parametrization of all the pp data [2£| and scale it for collision of lead nuclei. We account for the 
intrinsic kx of the partons by multiplying this scaled pp contribution by a factor of 2. We see that the sum of the 
thermal and prompt contributions provides a satisfactory explanation of the single photon spectrum at large kx in 
the present model. 

We also note that while an inclusion of the pionic bremsstrahlung improves the description of the data at the lower 
kx marginally, the theoretical calculations are still well below the estimated (see later) experimental results. We can 
think of two possible reasons for this shortfall. One possibility is that there may be additional reactions, e.g., involving 
baryons, which may contribute at low kx- The other, and possibly more plausible, reason could be the neglect of the 
pionic chemical potential |33j in the present work, and which can arise arise to- wards the end of hadronic phase |34j . 
Let us elaborate on this aspect. Recall that the pion-chemical potential (fi % ) may reach a value of about 60 MeV 
in collision of lead nuclei at SPS energies)^ when T drops to 100 MeV. This will imply an enhancement in the 
rate of production of photons from the hadronic matter by approximately a factor of exp(// 7r /T) for every pion in 
the entrance channel, and also some increase due to Bose enhancement factor in the exit channel, if there is a pion 
there. These aspects are under investigation. We may add however that the inclusion of these corrections should not 
drastically alter the following results, as the hadronic phase contributions, already dominate the yield of photons at 
low kx, considerably. 

The calculations at RHIC and LHC energies are performed in an analogous manner, taking dN/dy at y — as 1260 
and 5625 respectively for central collisions of gold (at RHIC) and lead (at LHC) nuclei. This procedure has been 
discussed repeatedly in literature; the only difference we have in the present work is that we use the state of the art 
results for the rates, and a more complete equation of state for the hadronic matter. 



C. Photon intensity interferometry at SPS energies 

As a first step, let us investigate the differences in the patterns of the intensity interferometry of photons having high 
and low kx- Thus we perform calculations by choosing kix, ipi, an d yi of photons such that two of the three momentum 
differences, q ou t, <?sidc, and qi on g, vanish in turn. This theoretical construction helps us obtain the corresponding 
correlation functions. 

We show typical results for k\x — 0.15 GeV, 1.05 GeV, and 1.95 GeV in Fig.0] We see a very interesting feature in 
the outward correlation function: as the transverse momentum of the photon increases, we see a clear emergence of 
two sources, one which has a smaller correlation radius (see larger <7 ut) and the other which has a larger correlation 
radius. This aspect has remained a recurring theme in results for photon intensity interferometry |10| and has its 
origin in the emissions from quark and the hadronic phases respectively which have vastly differing source dimensions. 
It is interesting that this feature noted in early exploratory studies has survived the vast improvements in the rate 
calculations and the dynamics of evolution. The emergence of the two source structure at larger kx is facilitated by 
the decreasing contribution of the hadronic phase there. 

This feature is not expected for pion intensity correlations as they leave the system mostly at the time of freeze-out 
only. 

The side-ward correlations show a decreasing radial dimension as kx increases, as expected for a transversely 
expanding source, and look Gaussian in nature. The longitudinal correlation function also shows a similar behaviour, 
as far as the variation of the correlation radius is concerned, however it is definitely not Gaussian in form. Thus 
neither the function (Eq|3 describes these variations satisfactorily, nor does the procedure of taking Ri = l/q~i, where 
q~i denotes the momentum difference where C — 1 has dropped by a factor of 1/e compared to its value at qi = 0, can 
do a full justice to these distributions. 
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FIG. 4: (Colour On-line) The outward, side-ward, and longitudinal correlation function for direct photons produced in central 
collision of lead nuclei at CERN . 



In view of the above, we have numerically evaluated: 

R^^l/<qf>; i= out, side, long (17) 

where 

2 /rfgig?(C-l) f , 

< q * >~ fd qi (c-i) (18) 
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FIG. 5: The transverse momentum dependence of outward, side-ward and longitudinal radii for photons from central collision 
of lead nuclei at CERN . 



and plotted it as a function of fciT in Fig- HI 

The kx variation of the i? ut reveals a very rich structure. In order to fully appreciate the observations, we recall that 
the source of photons has a dependence ~ T" exp(—E/T) where v > 2, and E is the energy of the photons in the rest 
frame of the fluid. This is very distinct from the corresponding dependence for pions, which varies as ~ exp(—E/Tf), 
where TV is the freeze-out temperature and contributions accrue from locations along the freeze-out surface. The 
"reclining chair" behaviour seen in this variation is an outcome of competition between a high temperature (limited 
to small radii and early times) and a large transverse flow velocity (predominant at large radii and late times) in 
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FIG. 6: The one dimensional correlation function for the kinematic window used in WA98 experiment [T3| . assuming a fully 
source and emitting only single photons. Central collision of lead nuclei is considered. 

producing photons having moderate kr- Thus the large radial flow ensures that photons having a large fcy can also 
be produced at later times, thus enhancing the so-called duration of the source measured as i? ou t — -Rside- We shall 
see later that the the 'seat' of this "reclining chair" gets narrower for RHIC energies and vanishes at LHC energies, 
as the increasing radial flow rapidly cools the system. A hint of this behaviour is also present in the momentum 
dependence of the side- ward correlation radius, which decreases almost linearly with fcx- The momentum dependence 
of the longitudinal correlation radius is completely different from the expected for pions [3(j • We see that 

^long f° r photons is almost inversely proportional to the transverse momentum. We believe that these differences do 
arise due to the difference in the source function for photons and pions. 

D. One dimensional analysis and comparison to WA98 results 

Let us now look at our results for the one-dimensional analysis of the correlation function in terms of the invari- 
ant momentum difference gi nv corresponding to the transverse momentum, and rapidity window used in the WA98 
experiment |13(. Our results with all the kinematic cuts are shown in Fig.[|J] In order to simulate the probabilistic 
selection of photons, we first generated a sufficiently large number of photons according to the thermal distribution 
calculated by us earlier, in the transverse momentum window of kx e [0.10, 2.5] GeV, and (randomly) distributed them 
uniformly over the azimuthal directions and the rapidity window corresponding to the experiment. Next we sampled 
pairs so that their average transverse momentum Kt was in the appropriate window. The correlation function was 
then calculated using the expression Eq.Q. The results were then averaged by binning in q mv . 

We see that our results are described to a reasonable accuracy by the form: 

C = 1 + 0.5 a cxp [-qf nv R- m /2] (19) 

where i?i nv » 6.5 fm for 0.10 < k^ < 0.20 GeV and 5.3 fm for the transverse momentum window 0.20 < kx < 0.30 
GeV. In order to compare our results with the numbers quoted the WA98 experiment |l3j| we need to multiply their 
results with \[2 as the fits discussed in that work do not include the factor of 2 in the exponential used in the present 
work fEa llO|) . Thus the corresponding experimental results for the i?j nv obtained by the WA98 experiment are 8.34 
± 1.7 fm and 8.63 ± 2.0 fm, respectively. While our predictions are in reasonable agreement with the "experimental" 
findings (within the errors), the experimental results are on the larger side. Before passing a judgment on our results, 
recall that that the values quoted by the WA98 experiment are obtained by assigning a values of 7.85 and 7.25 fm for 
the i?side, 8.4 and 7.7 fm for the i?i on g and 8.5 fm for the i? ut based on the values obtained for pion interferometry [3l| . 

Our results for the i? s idc are 4.93 and 4.85 fm, for the i?i on g are 3.7 and 2.5 fm, and for the i? ut are 9.97 and 9.39 
fm respectively at kr =0.15 and 0.25 GeV respectively (see FigEJ). Thus we feel that the results for the i?i nv as well 
as the results for the single photons at low kr quoted by the WA98 experiment (based on the value of the parameter 
a in the Eq. i|19l above) , are model dependent and uncertain to the extent that they use correlation radii determined 
from pion interferometry. We also see that at least the side-ward and the longitudinal radii for photon correlations 
are much smaller than the corresponding values for the pion interferometry. 
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FIG. 7: (Colour On-line) The outward, side-ward, and longitudinal correlation function for thermal photons produced in central 
collision of gold nuclei at BNL RHIC. 



A full three dimensional determination of the correlation function will go a long way in getting reliable results and 
constraining the theoretical calculations discussed here. 
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FIG. 8: The transverse momentum dependence of outward, side-ward and longitudinal radii for thermal from central collision 
of gold nuclei at BNL RHIC. 



E. Results for RHIC and LHC energies 

We now present results for correlation functions for thermal photons at RHIC and LHC energies (see Figsl7 HTT?|) . 

Looking at Fig. for central collisions of gold nuclei at RHIC energies, we note that the two-source aspect in 
the out-ward correlation function becomes more clear, due to increased contributions from the (smaller but stronger) 
quark-matter. Of-course, it is known |17| that we shall have a large contribution of pre-equilibrium photons at larger 
kx and this trend should continue at LHC energies. This will further enhance the contribution of the smaller-sized 
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FIG. 9: (Colour On-line) The outward, side- ward, and longitudinal correlation function for thermal photons produced in central 
collision of lead nuclei at CERN LHC. 



source, and the structure seen here, will disappear. 

The transverse momentum dependence of the source sizes( Fig. |SJ, is similar in nature to what we see at SPS 
energies, though the 'seat' of the 'reclining chair' seen in the out-ward size gets narrower, as one would expect for a 
more rapidly expanding source. The longitudinal correlation length is again seen to decrease rapidly with the increase 
in the transverse momentum. 

These trends continue at LHC energies (Fig. |5J) and the transverse momentum dependence of the source sizes 
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FIG. 10: The transverse momentum dependence of outward, side-ward and longitudinal radii for thermal from central collision 
of lead nuclei at CERN LHC. 



(Fig. I10[l becomes very pronounced. The so-called 'seat' in the outward correlation vanishes completely, and all the 
radii decrease roughly as 1/fc^, with a sw 0.27 for the outward correlation, about 0.14 for the side-ward correlation 
and about 0.94 for the longitudinal correlation. We note that at none of the energies, the longitudinal correlation 
function resembles a Gaussian, as normally assumed in parametrization. This is perhaps related to the emission of the 
photons from all the points in a longitudinally expanding system. In fact, similar shape is seen for the lon gitu dinal 
correlation function for pions 'escaping' from the system [9j and for pions produced from partonic cascades |36(. 
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FIG. 11: The difference of out-ward and side- ward correlation radii for thermal photons at CERN SPS, BNL RHIC, and CERN 
LHC. 

F. The life-time of the source 

The difference of the out-ward and side-ward correlation radii is often associated with the life-time of the source. 
In order to see the evolution of this parameter with the transverse momentum of the thermal photons we plot their 
difference as well as their ratios in Fig. 1111 for the three systems studied in the present work. We see the 'reclining 
chair' behaviour for the SPS energies, and a more rapid decrease as we go from RHIC to LHC energies. Several 
interesting results emerge. 

Firstly the difference (i? ou t — -Rsidc) is nearly unaltered as we go from SPS to RHIC for photons having transverse 
momenta less than 0.5 GeV. We have already noted that the 'seat' of variation is caused by the competition between 
cooling due to expansion and the blue-shift of the momenta due to the transverse flow. The more rapid decrease of 
the difference arises due to a more rapid cooling and an increased flow as the initial temperature rises. The same 
trend is seen for the ratios. An experimental confirmation or a refutation of theses predictions will be very valuable. 

III. DISCUSSIONS AND SUMMARY 

Before summarizing, it may be of use to discuss some aspects which we have overlooked in the present work. 

It should be of interest to examine the role of the reactions of the type q + q ^ "f + "f or q + g^q + "f + "f, for 
the intensity interferometry of photons. Firstly, the quark matter contribution to single photons is already marginal 
for low transverse momenta, which is the subject of our focus here. Then we note such photon pairs ( "diphotons" ) 
will be produced at 'a point' and not at different locations like a pair which comes from, say, q + q — > g + 7 and 
q + g — > g + 7 + 7. More-over their transverse momentum difference qx ~ 2kx, where kr is the transverse momentum 
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of either photon. For calculations which we can treat perturbatively, this difference is at least an order of magnitude 
larger than the values of qx where the correlation function has a significant value. Of-course one could include them 
as a source of single photons. However they are suppressed by a factor a s /a compared to the q + q g + J, say. 

Let us also discuss the role of direct (QCD) photons to the kx region under consideration here and for which we 
have experimental data at SPS energies from the WA98 experiment. In addition to the results given here, calculations 
have been reported by Dumitru et al. |37| . where large values for intrinsic momenta of partons are used to exhaust 
the WA98 data beyond kx ~ 3 GeV. The calculations exhaust only about 50% of the single photon yield at 2 GeV (as 
in the present work) and less than 20% of the yield at 1 GeV. Turbide et al. 20] have extrapolated the pQCD results 
down to zero transverse momenta (!) and find that the thermal photon yield is up to 4 orders of magnitude larger 
than the direct photon yield at lower transverse momenta. Parton Cascade Model calculations for SPS energies [14J 
show the dominance of quark matter contributions beyond 3 GeV (including photons which could be termed thermal, 
from the quark-matter) but a much smaller production from them at lower kx ■ These considerations should convince 
us that indeed thermal photons dominate in the region of transverse momenta considered here. 

We may also recall here the results of two works, which address similar issues. Peressounko ^lj has used hydrody- 
namics calculations with rates for the production of photons from quark matter and hadronic matter, along with the 
photons from the decay of pions to get A as well as the correlation radii, for direct photons. While the i? ou t and the 
i?i ong vs. kx behaviour estimated by him are similar to ours, his work shows a very unusual result: the i? s ide actually 
increases, and that too substantially, with increase in kx- Considering that photons having large transverse momenta 
are produced very early in the collision when the transverse expansion has not yet set-in, this is hard to understand. 
Similarly, the decrease in the values of the correlation radii, as kx decreases below 0.5 GeV at SPS energies is also 
difficult to fathom. Other differences could be attributed to the more complete rates for photon production used in 
the present work. 

The work of Renk ^l| is interesting in that it uses a fire-ball description to model the collision at SPS and RHIC 
energies. The parameters of the fireball are adjusted to reproduce the pion spectra as well as the HBT radii for pions. 
The differences seen in the present work then arise due to the explicit appearance of a mixed phase here, during which 
the speed of sound is zero and acceleration of the expansion is stalled. Results of an investigation using different 
equations of state, including one which does not admit a mixed phase, will be published shortly. 

In brief, the side-ward, out-ward, and longitudinal correlation functions for intensity intcrfcromctry of thermal 
photons at SPS, RHIC, and LHC energies are calculated. The transverse momentum dependence of these radii 
are very different from the corresponding results for pions, which are expected to decrease as 1/ \jmx for all the 
components. The longitudinal correlations for the three energies are quite similar and may be indicative of boost- 
invariance of the flow assumed in the work. The ratio i? ou t /-Rsido at LHC decreases rapidly and approaches unity, due 
to increase in i? s ide due to expansion and decrease in i? ou t due to rapid cooling. 

As the results are free from distortions due to final state interactions and uncertainties about the production- 
vertexes, or the production mechanism of the particles under investigation, a confirmation or refutation of the findings 
will be very valuable. 

It is also pointed out that the one dimensional analysis in terms of the variable Qinv may have only a limited use. 
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